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We continue to introduce bi-galileon theory, the generalisation of the single galileon model in- 
troduced by Nicolis et al. The theory contains two coupled scalar fields and is described by a 
Lagrangian that is invariant under Galilean shifts in those fields. This paper is the second of two, 
and focuses on the phenomenology of the theory. We are particularly interesting in models that 
admit solutions that are asymptotically self accelerating or asymptotically self tuning. In contrast 
to the single galileon theories, we find examples of self accelerating models that are simultaneously 
free from ghosts, tachyons and tadpoles, able to pass solar system constraints through Vainshtein 
screening, and do not suffer from problems with superluminality, Cerenkov emission or strong cou- 
pling. We also find self tuning models and discuss how Weinberg's no go theorem is evaded by 
breaking Poincare invariance in the scalar sector. Whereas the galileon description is valid all the 
way down to solar system scales for the self-accelerating models, unfortunately the same cannot be 
said for self tuning models owing to the scalars backreacting strongly on to the geometry. 
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I. INTRODUCTION 



Galileon theory was developed by Nicolis et al [l|, in order to facilitate a model independent study of 
certain infra-red modifications of gravity. They considered a class of scalar tensor theories of gravity, where 
all modifications of General Relativity are encoded in the Lagrangian for a single scalar field propagating 
in Minkowski space. The scalar field Lagrangian C(tt, dir, ddn) is invariant under a Galilean symmetry, 
7r — > 7r + b^x^ + c. The inspiration for the galileon description comes from co-dimension one brane world 
■^j- . models exhibiting infra-red modifications of gravity 0-111 • In these models, gravity on the brane is mediated 

by the exchange of the graviton and an additional scalar, often corresponding to the strongly coupled brane 
bending mode [6]. This strong coupling allows us to take a non-trivial limit in which the graviton and the 
scalar decouple in the AD effective theory Q. The scalar sector contains higher order self interactions and 
is Galilean invariant, a remnant of Poincare invariance in the original bulk spacetime. Many features of 
the original brane models such as self-acceleration Q , instabilities [9| and Vainshtein effects [ljj, [H| can be 
studied at the level of the corresponding galileon theory 0, [l2l - [l4j . 

In our companion paper [151. we introduced bi-galileon theory. This extends Nicolis et aFs model to two 
coupled galileon fields (see |la - [ig | for further extensions). Bi-galileon theory has particular relevance to 
co-dimension two brane world models exhibiting infra-red modifications of gravity [20l- [28| . Indeed, in [l5j . 
we showed that the boundary effective field theory for the cascading cosmology model [2CH22| corresponds 
to a bi-galileon theory in the decoupling limit. In an orthogonal paper fl8j . we considered a multi-galileon 
extension with internal symmetries, using the higher order interactions to evade Derrick's theorem and 
stabilise soliton solutions. However, in this paper, we shall return to the galileon as a means of modifying 
gravity. 

The most general bi-galileon theory |15l - il7| corresponds to a Lagrangian for two coupled scalar fields, n and 
£, propagating on Minkowski space. The Lagrangian is bi-galilean invariant, in that it remains unchanged 
under the following transformation it — s- tt + b^x^ + c , £ — > £ + b^x^ + c . By coupling one of the scalars 
directly to matter, through the trace of the energy-momentum tensor, we can interpret this as a modified 
theory of gravity mediated by the usual graviton plus two additional scalar fields. As in [lj, we neglect 
any direct mixing between the graviton and the scalars, although we do include mixing between the scalars. 
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This is consistent with the decoupling limit of the boundary effective theory for the cascading cosmology 
model. Indeed, we argued that this should also apply to any co-dimension two braneworld model exhibiting 
infra-red modifications of gravity. 

Although our previous paper [TEj was largely devoted to the formulation of our theory, we did initiate a 
study of maximally symmetric vacua, establishing elegant geometric techniques for assessing their stability. 
We did not, however, look at the phenomenological properties of these vacua. This is the subject of this 
paper. We will be particularly interested in two types of vacua: self-accelerating vacua and self tuning 
vacua. Roughly speaking, a self accelerating vacuum is one which undergoes cosmic acceleration even in the 
absence of any vacuum energy. There has been plenty of interest in these scenarios recently , since they 
represent a gravitational alternative to dark energy [231. We should be clear that self-acceleration does not 
go as far as solving the cosmological constant problem. One has to assume the existence of some unknown 
mechanism that sets the cosmological constant to zero, and then argue that the observed acceleration [3(| 
can be explained by new gravitational physics kicking in on cosmological scales. 

Self-accelerating vacua are often plagued by ghost-like instabilities ||. One of the aims of the original 
galileon paper [l[ was to see if self- acceleration is possible without ghosts. For a single galileon, ghost-free 
self acceleration is possible although to avoid overly restricting the domain of validity of the theory one must 
include a tadpole term which essentially renormalises the vacuum energy to non-zero values [if. In section 
IIIII we will study self- accelerating vacua in our bi-galileon theory, demonstrating the fact that ghost- free self 
acceleration is possible without any of the additional problems found in the single galileon theory. 

A self tuning vacuum is one which is insensitive to the vacuum energy. Such a vacuum is Minkowski 
even in the presence of a large vacuum energy. Self tuning mechanisms should, in principle, solve the old 
cosmological constant problem by dynamically tuning the vacuum curvature to zero whatever the vacuum 
energy. Whilst Weinberg's no-go theorem makes self tuning impossible for constant field configurations, 
it does not rule out more general scenarios (3l| . Indeed, co-dimension two braneworld models offer some 
hope for developing a successful self-tuning mechanism |20l - [27j . The reason is that adding vacuum energy 
to a co-dimension two brane merely alters the bulk deficit angle and not the brane geometry. Difficulties 
arise when one tries to study nontrivial branes geometries as this can sometimes introduce problems with 
singularities [24[ and perturbative ghosts [HJ . Furthermore, on a technical level, going beyond the vacua, 
even perturbatively, can be quite challenging in co-dimension two models. 

Given its association with co-dimension two braneworld models, it is no surprise that our general bi- 
galileon theory admits self tuning solutions. Indeed, as we will see in section ITVl the corresponding vacua can 
sometimes be stable. For both self tuning and self accelerating solutions, the study of non-trivial excitations 
is not too demanding. We will be particularly interested in spherically symmetric excitations so we can 
compare our results to observations within the solar system. Typically one would expect the additional 
scalar fields to ruin any agreement with solar system tests of gravity, mainly through the troublesome vDVZ 
discontinuity (33[. For both self acceleration and self tuning, we require the scalars to have an order one effect 
at cosmological scales and need some mechanism to suppress this on solar system scales. To achieve this we 
need to appeal to non-linear effects. Two possible mechanisms spring to mind: the chameleon mechanism 
[Hj], which relies on non- linear coupling to matter, and the Vainshtein mechanism [l(| HH, which relies on 
non-linear self interactions. Since we have assumed that there is a linear coupling to matter we must rely on 
the Vainshtein mechanism and scalar-scalar interactions. As we will see in section ITOI for self acceleration, 
this will help screen the scalars at short distances, so that our theory passes solar system constraints. In 
contrast, for self-tuning, we will see in section ITVl that one cannot "self tune away" a large vacuum energy 
and still hope to pass solar system constraints, at least at the level of the galileon description. This is because 
the backreaction of the scalars onto the geometry is too large and so the galileon description breaks down. 

Given a solution, our priority has been to establish whether it is perturbatively stable, and if so, whether 
it admits spherically symmetric solutions that can pass solar system tests of gravity through Vainshtein 
screening. For self-acceleration, this was in fact possible, but there are, of course, other things to consider. For 
example, what is the cut-off for our theory? This is particularly relevant because the Vainshtein mechanism 
is closely linked to strong coupling effects [11]. In DGP gravity, for example, fluctuations around the 
Minkowski vacuum become strongly coupled at length scales below around 1000 km, corresponding to a 
rather low momentum cut-off in the effective AD theory Q. However, one can argue that we do not live in 
the vacuum, so this is not necessarily a problem Q • The relevant cut-off comes from considering fluctuations 
about the non-trivial spherically symmetric solution arising in the solar system. This can sometimes be 
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much higher, as is the case for the asymptotically flat DGP solution [7(. Of course, we need to check this 
explicitly on a case by case basis. As the background solution changes with scale, so the strong coupling 
scale will run. There exists a critical radius at which quantum effects start to dominate and one cannot trust 
the classical background. We will require this scale to be shorter that the Schwarzschild radius of the Sun, 
below which we don't expect to be able to trust the galileon description anyway. 

Fluctuations about spherically symmetric solutions can also suffer from both superluminal and extreme 
subluminal propagation of fields. The latter can cause a large amount of Cerenkov emission, spoiling the 
quasi-static approximation of the solution [l]. Both of these effects were found to be problematic for the 
single galileon theory, but can be avoided simultaneously in the bi-galileon theory. Indeed, despite placing 
quite a few constraints on our theory, we will find that in some cases we tick all the boxes. This is in 
contrast to the case of a single galileon for which Cerenkov emission and a low scale of strong coupling 
can only be avoided in a ghost-free theory by introducing a tadpole. For self-acceleration, such a tadpole 
is theoretically undesirable as all it really does is renormalise the vacuum energy to a non-zero value. To 
further avoid superluminality problems, one is forced to eliminate all interactions making it impossible to 
satisfy solar system constraints. None of this is an issue in bi-galileon gravity - even without a tadpole, all 
of the would-be pathologies can be simultaneously avoided, including superluminality whilst retaining some 
higher order interactions and Vainshtein screening. This is really the main result of this paper: bi-galileon 
gravity can give rise to ghost free self acceleration without introducing any of the other pathologies that 
plagued the single galileon theory. This suggests that bi-galileon inspired cosmological models are worth 
developing further as a viable alternative to dark energy. 



II. THE BI-GALILEON MODEL 



Let us begin by reviewing the main results from our companion paper 15|, with some additional niceties. 
We considered a modified theory of gravity with two additional scalar fields. We treat the theory as an 
effective theory in Minkowski space with the following field content: a single massless graviton, h pv , and 
two scalar galileons, ir and £. Aside from the coupling to the energy momentum tensor, T pu , we neglect any 
interactions involving the graviton. Since we assume a linear coupling between the scalars and matter, we 
can assume, without further loss of generality, that just one of the scalars, 7T, say, couples to the trace of 
the energy- momentum tensor, T = r/^T^. If the n- matter coupling is extremely weak, there is very little 
deviation from General Relativity, rendering the theory uninteresting on cosmological scales. In any event, 
the general theory is given by 
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where ^f n are the matter fields, £ = 
galileon Lagrangian is 



with 
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. . is the linearised Einstein tensor, and the 
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We see from the matter coupling that the physical metric is given by gp V = rj^ u + h^, where h pi/ = 
hy, v + Given a source T^ v , gives the usual perturbative GR solution, and so 2-^7^ gives the 

modified gravity correction. The field equations for the scalars are 
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where a TO , n and b m ^ n are given by 

a m ,n = (m + l)(a m ,n + /3 TO +l,ra-l), 6 m ,„ = (n + l)(/3 m ,n + Om-l,n+l) • 

and satisfy 

na m -i, n = mb m , n -i' 



(5) 



(6) 



In [Ta |. we initiated a study of maximally symmetric vacua, where the background stress energy tensor 
can have a non-zero vacuum energy = —a 5£. For de Sitter space with Hubble radius H , we can 
expand the metric about the here (x = 0) and now (i = 0), so that for \x\ <C H~ x and |t| <C if -1 we have 
/i^„ = —^H 2 x a x a r]^. Motivated by this we took our background fields to be 



7r(x) — -rk^Xfi 



£(x) = --k^x^x 11 



(7) 



where fc^ = H 2 — Hq R is the difference between the Hubble parameters calculated in the modified gravity 
theory and in GR, and fcj is a constant. Now a particularly interesting quantity is the action polynomial, 
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which is closely related to the on-shell action for the background galileon fields, 
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0^m+ra<4 



+ 7f(-4(j) 



c? a; x^x^ 1 I L(kn, k(=) , (9) 



As we showed in [15] , a stable ghost- free vacuum is one that corresponds to a local minimum of the action 
polynomial. In other words, we require J^- — J^- = 0, for a solution. Furthermore, for the solution to be 

ghost free we require the Hessian, Hess(L)ij — g ® g k , to have non-negative eigenvalues. 

The stability conditions for the vacuum were found by considering perturbations on the background 
solution. We found that by shifting the fields tt —> ft + tt, £ —> £ + £, we get the same form for the field 
equations, but with different coefficients, 



(10) 



l<m+n<4 



As the vacuum is already a solution, there is no contribution from m = n = 0. For 1 ^ m + n 4, these 
new coefficients are related to the original ones by the linear map 



l m,n ~ / ! M mtn 



(11) 



where 



lvl m,n — I p. 



i-\-j—m — n 



(4 
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= z!/m!(i — m)! 



(12) 



For TV < 0, we extend the definition of "factorial" using the Gamma function, Nl = T(N + 1), and recall 
that 1/T(N + 1) = when N is a negative integer. It immediately follows that we only get contributions 
from i + j ^ 4. 

Given that the field equations have exactly the same form in the perturbed theory, it is clear that the 
effective action describing perturbations also has the same form, 
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where the coefficients in the equations of motion are related to the coefficients in the effective action via the 
standard relations 



a 'rn,n — ( m + l)( a m,n + Pm+l,n-l)t Kn,n ~ ( n + l)(An,n + a ro-l,n+l) ) 



(14) 



Note that we have also made a shift in the graviton field — > grou When we study spherically 

symmetric solution in section Till A[ we emphasise that we are working with the perturbed theory given by 
Eq. m 

Given an action polynomial for the background, we can reconstruct the corresponding action for the theory 
by computing the following coefficients: 



= (TO + l)(a m ,n + /Wl,n-l) = -(-2) 
(n + 1)(/? TO ,„ + Om-l.n+l) = = -(-2) 



,„(4-m-n)! d m+n+1 L 
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where m, n. = 0, 1, 2, . . ., and we recall that we define a_i, n = /? TO .-i = 0. Note that since tt£ m -i,n — ££m,n—i 
is a total derivative for n,m > 1, we are free to set, say, /3i )W = — (3%^ = ... = 0, without loss of 
generality. However, it is clear from Eq. fTS] that this choice subsequently fixes a m ^ n uniquely. 

The action polynomial can also be used to quickly calculate the equations of motion for fluctuations on a 
solution, given by Eqs. 1101 Given the action polynomial, we can compute the coefficients in the perturbation 
equations directly, using the following formulae, valid for 1 ^ m + n ^ 4 



\m+n 



-(-2) T 



(4-TO-n)! d m+n+1 L 
4!m!n! dk^ +1 dk^ 

, n (4-m-n)! d m+n+1 L 



A\m\n\ 



dk™dk™ +1 



(17) 
(18) 



We conclude our review with a word or two about backreaction of the scalars on to the geometry. 
In [l5j], we argued that this backreaction could be neglected provided T s ^ Iar [rj; n, £] <C M^Sh^, where 

^scaiarfe; = 77=^7 / d A Xy/=g £ sca iar and £ sca iar[ff; 7T, £] is constructed out of the covariant completion 



~ scalar 

7r , as described in [TBI ]. On maximally symmetric backgrounds, this condition 
schematically corresponds to 
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III. SELF ACCELERATION 



We begin our analysis with self accelerating vacua. A self accelerating vacuum is one that accelerates 
even in the absence of any sources for the physical fields, h^ v and tt. There is some ambiguity as to what is 
actually meant by this if tadpole terms are present. The point is that at the level of the graviton equations 
of motion, the source corresponds to the vacuum energy, a. However at the level of the scalar equations of 
motion, we note that the tadpole term, J d 4 x aQ^n, has the effect of renormalising the vacuum energy seen 
by the tt field, a — >• a + ao,o- To avoid "cheating", we set the bare vacuum energy a = 0, and require the 
7r-tadpole term to vanish, ao,o = 0. This is in line with our assumptions at the end of the previous section. 
It also guarantees that Minkowksi space is a solution for the physical metric since the field equations can be 
solved by = 0,ir = 0. Note that this Minkowski solution need not be stable. Indeed, our interest is in 
stable de Sitter solutions. Given the constraints a — 0, ao.o — 0, any de Sitter solutions are necessarily self 
accelerating. 
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Now, since a = 0, the corresponding GR solution is always Minkowski, so for any maximally symmetric 
solution, we have = 0. It is non-trivial solutions for the scalar 7r that enable self-acceleration. Recall 
that Tt{x) = —jk- R x fl x ,J -, £(x) = — jk^x^x^ , so for self acceleration we require a solution with fc w ~ Hq, 
where H is the current Hubble scale. For such a solution to be stable, it must correspond to a minimum of 
the action polynomial. It is easy to build a suitable action polynomial as we will now demonstrate. 

Let us assume that a stable self-accelerating vacuum exists, with k^ — H^, k^ = CHq, where £ is real and 
can be either positive or negative. We can build the action polynomial about this solution by performing a 
Taylor expansion: 



m+n 



4! 



-{K-Htr^-cHty 



(4 - m - n)\ 

(20) 

Note the lower value of to + n in the sum above; by starting the sum at m + n = 1 we ensure that 
k n = Hq, k^ — (Hq is indeed a solution. One can straightforwardly check using Eqs. [T71 IT51 and [T4l that 
the coefficients are indeed the same a' m n and f3' m n defined in the previous section. To avoid a ghost on 
this solution, we require the Hessian of L at k^ = Hq, k^ = C-ffp, to have non- negative eigenvalues, or 
equivalently, 



<io 



>0, 



fox > 0, a' w P' Q1 > 



«oi + P'k 



(21) 



As this corresponds to self-acceleration, we have set a = 0, and further require the 7r-tadpole to vanish, 

4! / r/2 \ m + n 
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(4 - m - n)\ 



i-f) C n [(m + l)< n + TO^„C] (22) 



The right hand side of Eq. (2D is a quartic polynomial in £. We can certainly choose the coefficients so that 
their exists at least one real root, thereby guaranteeing self acceleration. Perhaps the simplest example of a 
theory with a stable self- accelerating vacuum has the solution £ = 0, so that 
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l<m<4 
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(23) 



We can regard this equation as fixing a' 4 , having already chosen suitable values for a[ , a' 2 , and a' 3 . 
Indeed, let us assume that 



"m,n TT 2(m+n-l) 



9 9 

M a' - V 
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where u m ^ n and v m ^ n are dimensionless numbers of order one, unless otherwise stated. The condition ([23 
for vanishing n tadpole now takes the simpler form 
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(24) 



which we think of as fixing 1/4,0, having already chosen wi,o, 1*2,0) and 1*3,0- Furthermore, in order to guarantee 
stability, recall that we must choose Ui,0j Vo,i and Mo,i + i>i,o so that Eq. 1211 holds. 

Note that the overall scale fi will typically control the strength of the linearised coupling to matter, whereas 
the Hubble scale, Ho, controls the higher order ddir, <9<9£ interaction terms. The action, which can be derived 
exactly from the action polynomial by computing the bare coefficients a m ^ n and /3„ 
schematic form 



will now take the 
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(25) 
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Let us now consider the issue of backreaction. In order to trust our galileon description, we require that 
the scalar fields do not backreact too much on the geometry. For the self-accelerating scenario we have 
just described, the condition p^|) amounts to x 2 p?H^ -C M^Hq. Since we restrict attention to subhorizon 
distances and sub Hubble times, it follows that /i < M p i, which is consistent with similar conclusions drawn 
in the case of the single galileon field [1]. 



A. Spherically symmetric solutions and the Vainshtein effect 

One of the features of self accelerating vacua is that they require 0(1) modifications of General Relativity 
on cosmological scales. If such large deviations from GR still occur at much shorter scales, within the solar 
system, then it is clear that we will not be able to pass local gravity tests [36| . We need a mechanism to 
suppress the scalars in the solar system. We know of two: the chameleon mechanism 34], which makes 
use of a non-linear coupling to matter, and the Vainshtein mechanism which makes use of non-linear self 
interactions [35|. Since we have assumed a linear coupling to matter, we are forced to appeal to the latter 
to help screen the scalar degrees of freedom at short scales. 

Although the Vainshtein mechanism is intimately related to strong coupling in modified gravity theories 
[TlT |. it is a completely classical effect that is still not fully understood. The basic idea is that linearised 
perturbation theory around a heavy source breaks down at some "Vainshtein" scale, ry, below which non- 
linearities become important, helping to suppress fluctuations in the scalar field relative to the graviton 
fluctuations. Typically a heavy source like the Sun is treated as a point particle with a definite Vainshtein 
radius. This is, of course, an over simplification. The Sun is an extended object, made up of many point 
particles each with their own Vainshtein radii. Is the Vainshtein radius of an extended body the same as the 
Vainshtein radius of a point particle with the same mass located at the centre of mass? The answer is not 
known, although given the role of non-linearities in the Vainshtein mechanism, one mig ht expect the answer 
to be negative. A detailed study is beyond the scope of this paper, although see [37], [38| for some work along 
these lines. 

Having recognised some of the possible pit falls with this mechanism, we cautiously proceed in the usual 
way, treating the Sun as a point particle of mass M Q ~ 10 39 M p /. We will be interested in spherically 
symmetric profiles for the scalar field that asymptote to the self- accelerating vacuum, looking to establish 
the scale at which the linearised theory breaks down, and checking to see if the scalars are indeed screened 
below that scale. To this end, we will consider spherically symmetric fluctuations about a generic vacuum 
solution, which we will ultimately take to be the self-accelerating solution just described. Consistent with 
our model of the Sun, these fluctuations will be due to a point source of mass M located at the origin. This 
means we have 

7r = 7f(x) +7r s (r), £ = £(x) +£ s (r) 

where the background fields, tt £ are given by Eq. [7J and the fluctuations ir s (r), £ s (r) satisfy the field 
equations ([TP]) with a source 6T{f — MS 3 (r) diag(— 1,0). By explicit calculation, we have 

1 d , o , , 
r z ar 

P 2 d a 2d 
r z ar r z ar 

£3,0 = ir f s ) £2,1 = — -r{r-K s Q 

r z ar r ar 

£4,0 = £3,1 = £2,2 = 0, (26) 

and for m < n we take £ m ,n — £n,m\n^- Note that ir' s — d r n s and £' s = <9 r £ s . Defining y = Tr' s /r and 
z = £g/r, the equations of motion become 



^±[r 3 A(y,z)\ = MS 3 (r) C>7) 
r z ar 

1 d 



^[r 3 B(y,z)} = 0, (28) 
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where 



am = /r+2(/ 2 Q +/ 3 a ), 



B{y,z) = ft + 2(f 2 + f 3 ) 

n 



(29) 
(30) 



By integrating (f2"T)) and 



In / j a i,n-iV z -i 

i=0 i=0 

over a sphere, we can recast the equations of motion as two algebraic equations 

M 



~ 47 rr3 
B(y,z) = 0. 



In the linearised regime, we have 



and so 1 



a i,o a o,i 
b[,o bii 



yiin 
~-l in 



( M/iirr 3 



M 



zu. 



M 



det[Hess(L)} 47rr 3 ' det[H ess{L)} Airr 3 

We can easily solve these differential equations to arrive at the linearised solution for the scalars 

4& o,i M Mn ,,_, 46' li0 M 



(r) = - 



det[ff ess(L)} 47rr' 



det[Hess(L)] 4irr 



(31) 
(32) 

(33) 
(34) 
(35) 



To discuss the Vainshtein effect in detail, we now focus our attention on the self accelerating scenario 
described previously. In general we have b' Q1 ~ b[ Q ~ [i 2 , and det[Hess(L)] ~ [i 4 , which means 



(U 



\t 



(u 



Assuming the linearised approximation is valid, how much does the tt mode 



cause deviations from GR? To answer this consider the graviton fluctuation, corresponding to the standard 



Newtonian potential, (h^) 



N 



M 1 



, and compare it to the linearised field tt. 



(lin) 



(lin) 
7Ts 
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pl 



> 1 



(36) 



where the inequality follows from the absence of backreaction onto the vacuum. On solar system scales, 
deviations from GR are constrained to within one part in 10 5 [36| . Since deviations from GR are at least of 
order one in the linearised theory, and one has to resort to the Vainshtein effect to pass solar system tests, 
as expected. 

The Vainshtein effect kicks in when the linearised theory breaks down. Taking the generic scales this 
occurs when yu n ~ zu n ~ Hq. It follows that the Vainshtein radius is given by 



r v 



M 
H n 



1/3 



Mpi 



2/3 



M 



1/3 



10- g M 

KM 



1/3 



(37) 



Ideally we would like the Vainshtein radius around the Sun to exceed the size of the solar system. The solar 
system extends out to the Oort cloud, which is of the order 10 16 m from the Sun. Given that the Sun has 



1 Recall from our previous paper 15] that 



\Hess{L). 
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mass Mq ~ 10 M p i it turns out we can write r 



10 



1/3 



10 16 m, and so the Vainshtein radius 



is at least three orders of magnitude larger than the solar system. 

It is not enough to prove that linearised theory breaks down below the Vainshtein radius, we also need to 
establish whether or not the Vainshtein mechanism actually takes place. In other words, does the 7r mode 
get screened at distances r < Tyl Generically, at short distances, the cubic terms in Eqs. l29l and l29l will 
dominate, and we have 



It follows that for r < ry, we have Ivr. 



r 




(38) 



where 



Now let us evaluate 



M p i 



in the solar system, up to the maximum distance r Q 



(39) 



m. It follows that within the solar system 



10" "Me 

pi u 



1/3 



16 



N 



2/3 



10 



(40) 



Since deviations from GR should not exceed more that one part in 10 5 in this region (36|, let us take /i ~ M p i, 
so that we safely pass tests of GR in the solar system without running into problems with backreaction on 
the vacuum. 

Of course we should not only be worrying about backreaction on the vacuum, but backreaction on the 
spherically symmetric solution itself. For distances greater than the Vainshtein scale, the linearised theory 
holds and so the question of backreaction is dominated by the background, where it is known not to be 
a problem when /i ~ M p i. Below the Vainshtein scale, things are a little more complicated. By taking 



7T + 7T S 



HI 



■ r 2 y, ddir ■ 
that schematically, M^Sh^ 



H 2 



-V, with y 
M/r 3 and 



Hq^, along with the analogous results for £, we can show 



scalar 



m+n+1 



It follows immediately that 



T: 



(r v H ) 2 E °(!) 

l<m+n<4 



r J 



m+n—A 



< 1 



(41) 



(42) 



for r < ry. Thus backreaction of the scalars is not an issue, even below the Vainshtein scale, and we can 
trust our solution all the way down to the Schwarschild radius. 

Let us finish this section with a word on strong coupling. Strong coupling is inevitably linked to the 
Vainshtein effect [ll| on a given background. Therefore, since the Vainshtein mechanism takes place, we 
expect quantum fluctuations about the background vacuum to become strong coupled at reasonably large 
scales. To elucidate this let us consider the effective action (|13p describing excitations about non-trivial 
vacua. Focussing on the scalar sector, we first diagonalise the kinetic term, by performing a linear map of 
0(1) on the scalars, (tt,£) — > (7r*,£*). Schematically, the scalar action now takes the form 



S scalar [^* i £,* ] ^ n] — J d 



y(^) 2 - y (^*) 2 + + 0(1%)5T 



/i 2 J2 (0(l)dn t dn,+0(l)dnM*+0(l)dTT,d^)(—^ ] 1 " 



Hi 



(43) 
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where Ai and A2 are the eigenvalues ol the Hessian, Hess(L). Both eigenvalues are positive in a ghost free 
theory, and, for self accelerating vacua, of order /1 2 . Now let us canonically normalise the scalar fields by 
defining 7r = yAiT* and £ = V^C*? so that 



Scoter = / d A x - i(avr) 2 - i(9|) 2 + ~(O(1)tt + 0(1)£>T 

2 (0(1)9^ + 0(1)^ + 0(1)^) 1 (^S 



n-l 

I (44) 



2^m+n^4 



Here we see explicitly how /n controls the strength of the scalar coupling to matter. Given that this is Planck- 
ian, the scalars couple to matter with gravitational strength. Furthermore, it is clear that the interactions 
become strongly coupled at a scale A ~ (fiHfi) 1 / 3 . For fj, ~ M p i, this means that tree level interactions are 
only valid down to distances of the order Aq 1 ~ 1000 km, below which loop effects are important. This might 
lead one to question the validity of our classical description. However, it is important to realise that these 
scales correspond to modes propagating on the vacuum. Our real concern lies with quantum fluctuations 
propagating on the non-trivial solution around the heavy classical source 0]. We will study fluctuations 
about the spherically symmetric solutions in the next section. 



B. Fluctuations about spherical symmetry 



In the analysis of the single galileon [l[, promising solutions were found to exist, exhibiting self- acceleration 
without ghosts and with a suitable Vainshtein mechanism occurring within the solar system. However, 
fluctuations about the corresponding spherically symmetric profiles revealed a number of insurmountable 
problems. For example, at large distances, radial modes were found to propagate superluminally, leading to 
worries about causality. In contrast, at shorter distances, angular modes were found to propagate extremely 
slowly, so much so that the earth's motion through the solar galileon field would result in excessive emission 
of Cerenkov radiation, raising serious doubts about the validity of the static approximation. A third problem 
pertained to the domain of validity of the classical theory set by the scale of strong coupling on the spherically 
symmetric background. Whilst it was hoped that strong coupling would occur at higher energies than in the 
corresponding vacuum, the opposite was generically true. None of these pathologies could be eliminated in 
an entirely acceptable manner. The related problems of Cerenkov radiation and low scale of strong coupling 
could be alleviated by suppressing quartic and quintic interactions. However, to avoid a ghost in this scenario 
one needed to introduce a tadpole, which is undesirable as it amounts to a renormalisation of the vacuum 
energy to non-zero values. To further avoid issues with superluminality, one would have to do away with all 
higher order interactions and so abandon all hope of Vainshtein screening in the solar system. In summary, 
it seemed that a fully consistent self-accelerating scenario could not be found. 

In this section we will show that all of these pathologies can be simultaneously avoided in the case of two 
galileons. To see this, let us consider small fluctuations about the spherically symmetric solutions discussed 
in the previous section, 



€ = £> + 4> 



(45) 



Although we will ultimately be interested in the scale of higher order interactions, for the moment we will 
focus on the leading order theory given by the quadratic Lagrangian 

! ' ' ' (40) 



where the angular derivative dn 



dtp, and dot products between angular derivatives are 



understood. Note that the fluctuations form a 2 component vector $ = 
matrices take the form 



K = 



u = 




Ufa/, 


v = 


v <p<p 





















and so the kinetic mixing 



(47) 
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These matrices can be calculated explicitly by, for example, working out the equations of motion coming 
from (|46[) . then comparing this with the linearised equations of motion of the system. By using (|45[) and 
y = TT' s /r, z = £' s /r, the results are most elegantly expressed using (f3T)|) as follows 



1C = J+ -d r J, V=U 



-d r U 



(48) 



where 



U = Ei + 2E 2 + 2£ 3 , J = E x + 3E 2 + 6S 3 + 6E 4 



(49) 



and 



8 f a f) f b 
u yJn u yJn 



(50) 



where and /„ are given by Eq. 



Ghosts, tachyons and superluminality 



Now the resulting equations of motion describe the linearised perturbation theory, and take following form 



- 2 $ + \d r [r 2 Ud r <!>] 



vol $ = 



(51) 



This is all we need to study a number of issues, including the speed of mode propagation, as well as possible 
instabilities arising from ghosts and tachyons. Indeed, we can read off the "no ghost condition" immediately. 
We simply require that K. has non- negative eigenvalues, or equivalently 



> 0, 



> 0, 



To address the other issues we need to derive the dispersion relations for the normal modes of the system. 
This is a little involved since these do not correspond to <f> and ip if the cross terms , , are 
non-zero. To derive the dispersion relations for the normal modes, we first assume that the background 
changes slowly with radius compared to the fluctuations, so that we can treat it as roughly constant. The 
equation of motion becomes — ICdf $ + lAdf^ + V8q$> ps 0, where the radial variation in the mixing matrices 
is being neglected. We now take Fourier transforms 



[K,^-U V 2 r - Vpl] $(w,p r ,pn) w 



(52) 



where p r and pci are the momenta along the radial and angular directions respectively. It follows that the 
dispersion relations are given implicitly by det \JCuj 2 — hip 2 — Vpfi] = 0. Now it is convenient to write p r = 
p cos q,pn — psinq, so that the dispersion relations for the two eigenmodes can be written as w 2 = c±(q)p 2 . 
The sound speeds along the q direction are given by 



4M - i 



TrM(q) ± \J (TrM(q)) 2 -AdctM{q) 



(53) 



where A4(q) —K X U cos 2 q+K, 1 V sin 2 q. We now impose the condition < c\(q) < 1, which should be valid 
for all q. The lower bound prevents tachyonic instability, whereas the upper bound prevents superluminal 



mode propagation. We can make these conditions more explicit, 
is equivalent to requiring that M has non- negative eigenvalues 2 . 



The "no tachyon" condition, cj_ > 
In contrast, the "no superluminality" 



2 In fact, the eigenvalues of M. are just c 2 ^. 
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condition, c±(q) < 1 is equivalent to requiring that M. — I has non-positive eigenvalues. Writing M. 

Mas 

these requirements correspond to 



"no tachyon" M u + > 0, 
"no super luminality" + < 2, 



(1 - M w )(l - M^) > M H Mw 



(54) 
(55) 



Large distance behaviour 

Let us now study these conditions more closely. We begin with an analysis of the behaviour at large 
distances, r ~S> ry. In the single galileon case, the radial modes become superluminal at this scale [l|. To 
see whether the same thing happens here we calculate the matrices perturbatively. First, we consider the 



equations of motion Eqs. l3"Tlandl3"2"lorder by order, using the ansatz y = y®+y 



>')_ 



and y 



- rK 



+z (nl) + 



where "1" labels the "leading order" contribution, "nl" labels "next to leading order" , and so on. Note that 
we suppress the evaluation of the expressions | r >ry except for the final results. We find 



l 



1 



(56) 



which lead to the relations d r y"> ~ —3y^/r, d r y^ nl ^ ~ —6y( n[S> /r, ... (and similar relations for z). We also 
expand £„ = En + En ' + ... (except for Si, which is constant), and it follows that 



-3£ 



(i) 



-6£ 



(nl) 



rd r T,f 



-6£ 



(0 



(57) 



Notice that the next to leading term £3 is of the same order of the leading term £3' . Plugging this into 
our formulae Eqs. |48l [49] and [50] we find 



JC ~ £1 — 3£ 



(nl) 



6£ 



U ~ £1 +2£ 



(i) 



2(£ 



(nl) 



V « £1 - £^ - 4(£ 



-,(n£) 



.(On 



(58) 



where we have ignored 0(-p>) terms. To compute M. we need knowledge of JC 1 , which we also calculate 
perturbatively 



K,- 1 = Ei 1 +6£f x (£ 



("0 
2 



sr 1 



(59) 



and so 



X = / + (3cos 2 q-l)£ 1 1 4° 



(6cos 2 g + 2)£ 1 l Ej 



(0 



(60) 



Now by studying the leading order contribution, it is immediately clear that, generically, the eigenvalues 
of M. — I will change sign at cosg = ±1/a/3, meaning that propagations along some directions will be 



superluminal. This can be avoided if we impose the condition £ 



(i) 



_ v (ni) 



r»rv 



v -i v (0 



= and require 



r>rv 



has negative eigenvalues. This guarantees that M — I has negative eigenvalues to leading 

order, ensuring that all modes propagate subluminally. We also see that there is never an issue with 
tachyonic instability at large distances since the sound speeds are always close to unity. Furthermore, as 
long as the vacuum is ghost free, then we have that JC has non-negative eigenvalues to leading order (£1 is 



non- negative definite), and therefore there are no ghosts and the condition that Ej *£ 



definite reduces to that £ 



(0 



r>ry 



is negative 



is negative definite. 
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Short distance behaviour 



We now turn our attention to shorter distances, below the Vainshtein radius, r -C ry ■ For the single 
galileon, short distance modes generically propagate very slowly along the angular direction. In fact, the 
propagation is so slow that one gets excessive emission of Cerenkov radiation as the earth moves through the 
solar galileon field, raising serious doubts as to the validity of our static approximation. This problem can 
be avoided by eliminating quartic and quintic interactions of the single galileon, although in the absence of 
a tadpole, this procedure renders the self-accelerating vacuum unstable to ghostly excitations. We will find 
that such difficulties can be avoided in the case of two galileons, even without having to introduce a tadpole. 

We first do the perturbative expansion at short distances: y = j/W + j/"') + y = + z^ + ... and 
£„ = + E^ + Y^k nV> + .... Note again that we suppress the evaluation of the expressions | r <r v except 
for the final results. By considering the equations of motion order by order, we find 

yQ\ Z W~l, i^O^nD^l, ... (61) 



r 



and 



zj 4 ~ -3, ^4 ,^3 ~ -j, ^4 ,^3 ,^2 >•■• \ VZi ) 

Similar to the large distance case, we can eliminate rd r E„ in favour of E„ and get 

K w 2£< n,) + 4E< n nl) + 2E< () + 4££ n,) + 2E<° , W « 2£< + 2E<"° + 2E 2 ° , V « 4"° + £<° (63) 

where we have neglected 0(1) terms. As we shall see in next subsection about the strong coupling issue, 
we can set E 4 = (no evaluation at r <C rv), i.e., set a' mA _ m = b' mA _ m = 0, to avoid the strong coupling 
problem. We shall assume this here and the three matrices reduce to 

/C« 2E^ } +4E3" + 2£ 2 °, 2E^ } +2E^"° + 2£ 2 °, VsEf' + sf 1 (64) 

Now, assuming det £3^ 7^ 0, we have /C" 1 « 5 [£3 and so ~ cos 2 q/. It is clear that the sound speed 
along the angular direction, co&q = 0, will be very small, and we will once again run into problems with 
emission of Cerenkov radiation due to the earth's motion through the slowly propagating galileon field. So 



we enforce the condition E 3 



(0 



= 0. Then we have 

r«rv 



M = (cos 2 q + 1)(4£ 3 "' } + 2E 2 () )- 1 (E^ } + £ 2 ^) + . . . (65) 



and we should also impose the condition that 
values. 



(4E^ ) +2E 2 ' ) )- 1 (E^ i) +E 2 °; 



has non-negative eigen- 

r<rv 



Strong coupling 

We saw in the previous section that quantum fluctuations on the self accelerating vacuum become strongly 
coupled below a length scale of around Aq" 1 <~ (/iiJ, 2 ) -1 / 3 <~ 1000km, beyond which one cannot trust the 
classical description. We are now ready to ask whether or not the corresponding momentum scale can be 
pushed higher for quantum fluctuations on the spherically symmetric solution. Generically this was not the 
case for a single galileon, and although the same is true here, we will find examples where there are no issues 
with strong coupling, or indeed any of the other pathologies we have discussed. For the single galileon, the 
scale of strong coupling can only be raised at the expense of introducing either a ghost or a tadpole, neither 
of which is theoretically desirable. Such difficulties are avoided in bi-galileon theory. 

Since strong coupling kicks in at short distances, we work in the short distance approximation. To 
estimate the strong coupling scales above the spherically symmetrical background, we need to work out the 
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interaction terms in the Lagrangian. While it is necessary to keep the background coordinates exact, ie, 
spherical coordinates, it suffices to treat the perturbed coordinates as cartesian and treat the coefficients as 
constants, as long as the strong coupling length scale is smaller the radius of the point in question. It is 
easier to work at the level of the equation of motion and then promote it to the Lagrangian. Neglecting the 
details of the contractions, and taking it' s /r = y to imply ddir s ~ y, along with the analogous results for £, 
the full Lagrangian is schematically given by 



+ «X E E {^r^\o m , + {m {^)\%' (66, 

2<m+n<4 2<i+i<4 \ ' \ U / \ U / 

where w ~ y ~ z, as given by Eq. [38J Now let us canonically normalise the scalar fields by performing the 
linear map <f> = iJZL (o(l)$ + 0(1)4>) , V = ^ \/¥ + °Q)$) > 



£ w «0(l)(^) 2 + 0(l)(^) a 

E E fe) A o (1 "" J) [0(1)^ + 0(1)^] (^) '(W)' (67) 



»+n-|(i+j)- 

(7k) ' 

2<m+n<4 2<i+j<4 

In (hopefully) obvious notation, these interactions become strongly coupled at a scale 



2(m + r»)-3(i+j)-l 

W 



6(l-i-3) 



A m +„,i+j ~ A ( ) (68) 



Since w 3> i?o a ^ short distances, it follows that the largest interaction comes from the term with i + j = 2, 
m + n = 4, becoming strongly coupled at 

A 4 , 2 - A (^J , (69) 

Thus strong coupling kicks in at an even lower scale than in the vacuum if such interactions are present. 
This is clearly undesirable, so we set all terms with m + n = 4 to vanish, a! m 4 _ m = b' m 4 _ m = 0. This does 
not come into conflict with any of our previous phenomenological constraints. Things are now much better 
behaved. The largest interaction comes from the term with i +j = 2, m + n = 3, becoming strongly coupled 
at a much higher energy scale 

/ \ 1/6 

A 3 , 2 ~ A ( , (70) 



H 



— p- 2 - J , from which we see that the strong coupling 

scale runs, and is given by 

A sc (r) ~ lO-WAV'r- 1 /* (71) 

Our classical description breaks down when rA sc (r) ~ 1, after which quantum corrections become important 
Q- This occurs at a critical radius r c ~ 10~ 13 ^ 5 A " 1 ^2 — 3 km, which is of the order the Schwarzschild 
radius of the Sun. This is perfectly acceptable from a phenomenological point of view as we don't expect 
the galileon description to be valid at such a low scale anyway! 
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An explicit example 

Having established the conditions, we point out that it is not difficult to engineer a bi-galileon model that 
achieves self acceleration without running into the various pathologies that hampered the single galileon 
case. Here we provide an explicit simple example that satisfies all the conditions: a\ = 3/z 2 , a' 1 = —^ 2 , 
&0,i = ¥ 2 , 4,o = V/Jtf, 4,1 = -13/i7#o> 4,2 = 24/x 2 /^, < 3 = -9^/H$, b[ l3 = -I8 f i 2 /H$ 
and 4 _ m = b' mi _ m — a' m 2 _ m = b' m ,2-m = 0- (To make contact with the overlying theory, we simply 
reconstruct the action polynomial from these parameters, using Eqs. [T41 and [20l) As the reader can check, 
there is no ghost in this model. Vainshtein mechanism is effective owing to the existence of higher order 
galileon terms. Since £4 = 0, we do not have the strong coupling problem. At small distances, since 

0, we have A4 ~ (cos 2 q + l)/4, and so the sound speed is always between 



^2 



and £ 3 ° 



I ^ c «L<rv ^ 5' ^ nus there is no problem of Cerenkov emission or superluminal propagation. At large 



distances, since £2 = and £3 is negative definite, the propagation speed is always around 1 and 

r>rv 

subluminal. Note that we have made use of the equations of motion. To connect to the dark energy problem, 
we can set Hq to be the Hubble scale and fi ~ M p i, in which case there is no problem of backreaction. 



IV. SELF TUNING 



Now let us turn our attention to self-tuning vacua. We define a self tuning vacuum as one that is Minkowski, 
for any value of the vacuum energy a. Self tuning mechanisms are designed to solve the "old" cosmological 
constant problem 3 . This arises because each of the matter fields, fy n , contribute to the overall vacuum 
energy density, er = X)n(P*n)- Each contribution is found by summing up the zero-point energies, 
of the normal modes, up to some cut-off A 3> m, where m is the particle mass. Setting H = c = 1, this 
generically gives [3l[ 

. , f A d 3 k 1 fZ - A 4 

so that the overall vacuum energy a ~ A 4 . We often assume that this cut off is Planckian giving rise to a 
huge vacuum energy of around 10 72 (GeV) 4 . We can reduce this slightly by allowing for supersymmetry, so 
that the cut-off corresponds to the supersymmetry breaking scale, but even then the vacuum energy is at 
least (TeV) 4 . In any event, in GR, such a huge vacuum energy would cause the vacuum to be highly curved, 
giving a de Sitter or anti de Sitter geometry with curvature scale ~ y/\a\/M p i. In terms of our graviton 

mode, we have h va J- = g§rx a x a r] tll/ . 

pi 

For self tuning to occur, we require the physical metric to be Minkowski, whatever the choice of a. In other 
words, we have h^ v = 0, and so 7f = — 12 Zr'i x a x a in order to cancel off the contribution from the graviton. 

It follows that we need — —j§p- to be a solution to the background field equations, j^- = J^- = 0, for 
any value of a. 

It is important to realise that the coefficients a m n , f3 m>n do not depend on a, as these coefficients define 
the theory and are independent of the source. Having said that, there is explicit a dependence in the action 
polynomial © 

N ^ / 41 

L(fc 7r , % a) = 4aK ~ V {u m .nK + P m , n h) -- Ta rrC fc e ■ ( 73 ) 

e -^ 1 \ 2 1 (4 — m — nil s 



3 See [39l for a nice review of the "old" cosmological constant problem. 
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and so different values of a give different action polynomials. For any given o, we require that the corre- 
sponding action polynomial has a minimum at k n = — , so that our self tuning is stable against ghost-like 

pi 

excitations. 

Now if the action polynomial L^k^jk^a) has an extremum at k n = — F.^ , then it is obvious that the 

pZ 

modified action polynomial 



L(k 7I , fej) = L(k 7T7 k^;a) — 6M 2 

pi J :jlv± pi 



3M3 9M 4 



(74) 



also has an extremum at this point. However, the interesting thing to note is that L is independent of cr. It 
follows that the modified action polynomial has a continuum of extrema, (k v , fef) = (C, /(C)), parametrized 
by C = — <r/3Mpj. Now, this is the case if and only if L has the form 

L(fc,r, fcj) = (&£ - /(fc 7r )) 2 C(fc 7r , %) + constant (75) 

where C(k v ,k^) is a cubic in fef, with fc^-dependent coefficients. As the constant is irrelevant, we might as 

well neglect it for simplicity. We can now think of L as a quintic polynomial in k(=, with a real double root. 

To establish whether or not we have ghosts, we need to compute the Hessian of the action polynomial on 
the solution. It is easy to check that 

Hess(L) = Hess(L) + ( 12 ° Q \ (76) 

It follows that at (k n , k^) = ((, /(C)), we have 

Hess(L)\ {u{c)) - ^ -2f(0C({J(0) 2C( C ,/(C)) J (7?) 

To avoid ghosts, we require the eigenvalues of this Hessian to be non-negative, so we simply need C(£, /(C)) ^ 
0. Therefore, the action polynomial for a theory admitting stable self-tuning vacua can always be written as 

5 

L(k„, % o) = 4ak„ + c 5 _„(fc 7r )(fc c - f{K)) n + QM^kl + constant (78) 

where ca(k n ) ^ for k w G ( — dgr, — ^frr ) • Such a self tuning theory is stable for a range of vacuum 

\ pi pi / 

energies a G (02,01). The function Lik-^^k^) ought to be a bivariate polynomial, up to fifth order in the 
variables k n and k(=. For generic functions /, c„, this will not be the case given the form for L in Eq. 1781 
Thus / and c n ought to be chosen appropriately. 

As a simple example, we can choose f(k n ) = M 4 /k 7T and c 3 (fc Tr ) = k^/fj. 2 , with other cs_„ vanishing. This 
gives rise to the self-tuning model: 

M 4 1 

£^ = 3M p V& - — ttOC + 5-5* [Dn((n0 2 - {d^f) - 2d ll d v nd»d v &3£ + 2d„d u 7rd 1/ d p ^d p d^] (79) 

One might check explicitly that the background equations of motion are solved by k^ = M A jk^ = 
—3M 4 M 2 l /a, irrespective of the value of vacuum energy a. Whilst we have a ghost on the trivial back- 
ground (k-n — = 0), this is not important since fluctuations on the self-tuning background are ghost free, 
as one can easily check. The backreaction of the self-tuning background onto the geometry will be negligible 
if the condition M 4 < /j,M p iHq is satisfied. Vainshtein screening is effective owing to the presence of the 
higher order galileon terms. Indeed, by performing perturbation analysis on the linear spherically symmetric 
solution, we can estimate that the Vainshtein radius of the Sun is about ry ~ (MqM p iM 4 / fia 2 ) 1 / 3 . Un- 
fortunately, as we will see in Section liV B[ one cannot self-tune a large vacuum energy without introducing 
problems with backreaction. 
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A. Evading Weinberg's no go theorem 

Let us divert our discussion and consider for a moment our self tuning solution in the context of Weinberg's 
no go theorem 31]. Weinberg argued, on very general grounds, that no dynamical adjustment mechanisms 
could be used to solve the cosmological constant problem. Let us briefly sketch his proof, adapted to the 
case at hand. Imagine a system of two scalar fields, 7Ti,7T2, non-minimally coupled to gravity, described by 
a general action 



(80) 



We assume that the matter fields, all lie in their ground state, and absorb their contribution to the 
vacuum energy into the potential for the scalar fields, V(tti, 7^). Let us consider a Poincare invariant solution 
to the field equations, with constant scalars, and "constant" metric, <? M „ — rj^. It follows that this solution 
satisfies the equilibrium condition 



dL 



dg^v 



£/,7r— constant 







g,7T— constant 



(81) 



For this solution to be "natural" we demand that the trace of the gravity equation is a linear combination 
of the scalar equations, 



9V 



dL 
dg^u 



^ OTTi 



(82) 



for all constant fields. This ensures that the trace of the gravity equation vanishes automatically by virtue 
of the scalar equations of motion 4 . Defining x — 7r i/2/i — 712/2/2 and <f> = 7ri/2/i + 112/2/2, we note that 
(for constant fields) the Lagrangian L is invariant under 



8(j> 



(83) 



It now follows that the Lagrangian must take the form L 
and so 



-g£(x, dXi dg^, dd . . .) where g^ v = e^g^, 



dL 



dguu 



g,7v— constant 



" n P-V T I 

~.y Js\g.7T= const ant 



(84) 



Applying Eq. 1811 we find i| g . w=C onstant = 0, which fine tunes the potential V(7Ti,7r2) to be vanishing in the 
Minkowski vacuum, ruling out a solution to the cosmological constant problem by dynamical adjustment of 
the fields. 

One might hope to promote the self tuning theories in the bi-galileon model to fully covariant theories that 
evade Weinberg's no go theorem. Although it is conceivable that self tuning is spoilt by covariantisation, one 
ought to recover it in the decoupling limit M p i —> 0. Assuming for the moment that it is not spoilt, it is clear 
that we evade Weinberg's no go theorem by breaking Poincare invariance. On the self tuning background, 
the physical metric is certainly Minkowski but the scalar fields are not all constant, rather 7r oc x^x^ . Of 
the Poincare symmetries, only translational invariance is broken, while Lorentz invariance is preserved. 



B. Spherically symmetric solutions and the breakdown of the galileon description 



Of course, we do not live in a vacuum, and so it is important to ask what happens when we introduce 
a heavy source into our system. Here we are interested in spherically symmetric excitations of the self 



4 Eq. 1821 is a sufficient but not a necessary condition for this to hold. 
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tuning vacua around the Sun. As in the self accelerating scenario, we need some mechanism to suppress 
modifications of GR at short scales. Here this corresponds to the break down of linearised theory through 
the Vainshtein mechanism. The equations that govern the excitations are the same as those given in section 
lIlTAl by Eqs. [23 [251 M and M 

However, it turns out that although one can engineer a Vainshtein effect at the level of these equations 
around self tuning vacua, one cannot do so without introducing a large amount of backreaction and destroying 
the galileon description altogether. To see this most efficiently, we will present an heuristic argument that 
illustrates the problem succinctly. Recall that the scalars will backreact heavily on to the geometry unless 

\TZ la r[r,^Z}\<<M 2 pl \£hn (85) 

where W v = W v + 2-K'rf- 1 ' is the physical metric perturbation and 

2 S 



9 5g^ v 



d xy/—g £ s 



-'scalar 



Here £ S caiar[<7; 7T, £] is constructed out of the covariant completion of — C^.t — SM^itd tt , as described 
in . The full set of galileon equations can be expressed as 

E *M = s^J d4x = - ? ^ TA1!/ ' E &>® = si J c ^ = °' EUW = m 2 ] (86) 

where tt — 7f+7r s (r), £ = £+£ s (r) correspond to the scalars evaluated on the spherically symmetric excitation 
about the self tuning vacuum, and h^ v ~ h^ c + h% u is the corresponding graviton. The energy-momentum 
tensor has two pieces: a large vacuum energy and the contribution from the Sun, 

T^ v = -orf v + 6T£ U (87) 

For self tuning vacua, the physical metric only contains a contribution from the spherically symmetric 
excitation hf" 1 — h% v = h^ v + 2-K s rf ,v . This is because the vacuum contribution to the physical metric 
vanishes on account of h^ c being cancelled by 2^1]^ . Furthermore, in the event of a successful Vainshtein 
mechanism, the graviton excitation should dominate the scalar at short distances, so we have W v « h% v . 

Since the equation for the graviton excitation is really just £ h^ v — , i| , it follows that at short distances, 

P i 

below the Vainshtein scale, 

M 2 vl \Eh^\ = \5Tg v \ (88) 

This fixes one half of the inequality 1851 governing backreaction. We now turn our attention to the other 
half, by first noting that, as a result of diffeomorphism invariance, 

d ll TZ l Jv^^)=^n + E 6 d^ (89) 

where 

From the galileon equations of motion (1551) . it follows that on shell: E^ = — t^JXq 1 ' — 6Mp ; D7r s and E^ = 0. 
In the Vainshtein region, linear contributions are subleading, and so one can neglect the linear 7r s term in 
E,,., and approximate it as E T ss —"H^ST^. Now plugging all of this into Eq. (S3 we find that 

Note that we have used the fact that ir ^> ir s . This is certainly true except at extremely short distances of 
the order \x\ < M p ij ^J\a\. For a large vacuum energy, this short distance cut-off is tiny, corresponding to 
around a millimetre for a TeV scale vacuum energy. 
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In terms of scale, our analysis suggests that 

|rr a laJ^]M^l|«T| (90) 

For the inequality (I55|) to hold, it is clear that we must have |7f| <C 1- However, recall that n — — 12 f^ x v x ^> 

which means that tt is large on solar system scales for large vacuum energy. For TeV scale vacuum energy, 
this suggests we have a big problem with backreaction at solar system scales. In contrast, for a ~ (meV) A 7 
the backreaction is small on subhorizon \x\ < H^ 1 , although there is little phenomenological motivation to 
self-tune such a small vacuum energy. 

Our conclusion is that one can self tune away the vacuum energy in the bi-galileon model; but if the vacuum 
energy is very large (e.g. TeV scale or larger) as predicted by current particle theories, it is impossible to 
do so without abandoning either the Vainshtein effect, or the galileon description, in the vicinity of the 
solar system. Of course, our arguments are suggestive rather than precise, so one might wish for a subtle 
resolution of this problem. We are, however, pessimistic. We have been unable to find a numerical example 
that does not behave in precisely the way suggested by our heuristic argument. 



V. DISCUSSION 

We have studied interesting phenomenological solutions to the bi-galileon model [l5[ • This can be under- 
stood as the decoupling limit of a gravity theory in which GR is modified by the addition of two scalar fields 
taking part in the gravitational interaction. This model is an extension of the original galileon model [l| to 
two galileon fields, and is expected to be particularly relevant to co-dimension 2 braneworlds [l5j . 

We have focused on two particular types of solution: asymptotically self-accelerating solutions and asymp- 
totically self tuning solutions. Let us first comment on self acceleration. In contrast to the single galileon 
case we have shown that one can find bi-galileon theories that do not contain any tadpoles, and admit self 
accelerating solutions that satisfy each of the following: 

• fluctuations about the vacuum do not contain a ghost 

• spherically symmetric galileon fields undergo Vainshtein screening in the solar system 

• fluctuations about the spherically symmetric galileons are never superluminal 

• fluctuations about the spherically symmetric galileons never lead to trouble with excessive emission 
Cerenkov radiation. 

• do not have an unacceptably low momentum scale for strong coupling, leading to the breakdown of 
the classical solution in the solar system due to large quantum fluctuations. 

• do not suffer from problems with backreaction, leading to the breakdown of the galileon description 
for either the vacuum solution or the spherically symmetric solution. 

It is, perhaps, remarkable that we can simultaneously achieve each of the above in a given model, in contrast 
to what could be achieved for the case of a single galileon. We believe this merits much more investigation 
into the bi-galileon model as an alternative to dark energy. 

Of course, the first step is to promote a good theory to a fully covariant one. There are two ways in 
which we might think about doing this. The first is to simply take our effective AD theory and perform a 
covariant completion, along the lines described in [13, S3] ■ Although, the Galilean invariance is broken, we 
expect the generic features of the galileon solution to be retained, at least up to corrections which are Planck 
suppressed. An alternative, more ambitious, approach would be to try to oxidise our theory and interpret it 
as a particular co-dimension two braneworld model with very desirable properties. To this end we note that 
the probe DBI brane description [4lJ was very recently generalised to the case of multi-galileons (l9j . 

One might reasonably ask how natural our "good" theories are? How stable are they against radiative 
corrections? Radiative corrections will typically come from two different sources: (i) galileon loops that will 
renormalise the coefficients in the action; and (ii) matter loops that can potentially introduce non-galilean 
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invariant terms as the coupling to matter breaks the Galilean symmetry. If we treat the galileon theory as 
an effective theory valid up to the strong coupling scale, A sc , we do not see any problem with naturalness 
arising from galileon loops. In contrast, matter loops are potentially more dangerous, as the effective theory 
for the matter sector is valid up to nearly a TeV. Of course, this is at the origin of the old cosmological 
constant problem, and its resolution is beyond the scope of this paper. 

This brings us nicely on to the other kind of solution studied in this paper: the self-tuning solution. By 
breaking Poincare invariance one can escape the clutches of Weinberg's no theorem |3lj . such that in the 
presence of a vacuum energy, the scalars adjust themselves accordingly, and eliminate the resulting curvature. 
The problems start when one requires self-tuning of a large vacuum as predicted by current particle theories 
and tries to study spherically symmetric solutions sourced by the Sun. Although it is possible to engineer 
Vainshtein screening at the level of the galileon description, one cannot do so without the galileon description 
itself breaking down due a large amount of backreaction of the scalars onto the geometry. We anticipate 
that this will make it extremely difficult to satisfy solar system constraints in a covariant completion of our 
self tuning galileon models. 

Intuitively, this is actually quite easy to understand. On cosmological scales, we are asking the galileon 
fields to do an awful lot of work in screening the large vacuum energy from the resulting curvature. Indeed, it 
requires a background galileon field, 7f| se if-tun 3> 1 at Hubble distances, in contrast to self-acceleration which 
has 7r| S eif-acc ~ 1 at Hubble distances. At the same time we are asking that the galileon fields do nothing on 
solar system scales, that they are screened by the Vainshtein mechanism and one is able to recover GR. It 
seems that such a dramatic change in the galileon behaviour between cosmological and solar system scales 
is impossible to achieve. 

As we stated in the introduction, the Vainshtein mechanism is not the only means of suppressing scalar 
fields in the solar system. If we are prepared to break Galilean invariance in the vacuum theory then 
we might consider an alternative mechanism whereby the scalars develop a large mass in the vicinity of 
heavy objects like the earth. Such mechanisms include the chameleon [34[ and the symmetron [42j], and 
we could even consider using them in tandem with the Vainshtein mechanism. Even so, as described in 
the previous paragraph, we are asking the scalar fields to change their behaviour dramatically, perhaps too 
dramatically. By making use of chameleons, we might also have to worry about possible violations of the 
Equivalence Principle [38|. In the event of an unsuccessful resolution of these issues, it would be worth 
asking if Weinberg's no go theorem [3l[ can be extended to allow the breaking of Poincare invariance, but 
taking into account phenomenological constraints coming from local gravity tests [36j . 
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